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1 Introduction 



One of the possible ways for building completely integrable quantum systems is the inverse 
method of separation of variables based on the use of the so-called multi-parameter spectral 
equations. The idea is to interpret these equations as the result of separating variables in a 
certain multi-dimensional completely integrable quantum system and to reconstruct the form of 
the latter by eliminating some of the spectral parameters considered as separtion constants (for 
details see e.g. refs. I], ||, 0, § ) • 

If one wants to obtain this way exactly solvable completely integrable systems, one should 
start with exactly solvable multi-parameter spectral equations. In this paper we present one pos- 
sible way of building such equations by considering as an example one-dimensional second-order 
differential equations. We show that in this case both the multi-parameter spectral equations 
and their solutions can be constructed from one and the same elementary building blocks called 
in the paper the p-functions which satisfy special functional relations which we called scalar 
triangle equations. The choice of this terminology is not accidental: The reader having some 
elementary acquaintance with the celebrated r-matrix approach (see e.g. |||]) can easily be con- 
vinced that there are many common features between the scalar triangle equations and ordinary 
matrix triangle equations^ There is also a deep relation between the solutions of these equa- 
tions: the /o- functions are natural analogues of the classical r-matrices. The general form of 
the substitution for solving the multi-parameter spectral equations is very similar to that of the 
famous Bethe Ansatz used for solving completely integrable quantum models. The numerical 
equations determining the solvability conditions have the same meaning as the corresponding 
Bethe Ansatz equations. 

In this paper we show that there are two essentially different classes of exactly solvable multi- 
parameter spectral equations which can be constructed from p-functions. We call them rational 
and irrational equations stressing the fact that in a certain " canonical" coordinate system their 
"potentials" can be expressed respectively via rational and irrational functions. The rational 
multi-parameter spectral equations are very well known in the literature. In 1987 Sklyanin 
obtained such equations Q as a result of a separation of variables in the so-called completely 
integrable Gaudin models ||T| (which are easily constructable in the framework of the r-matrix 
approach). Applying to the rational multi-parameter spectral equations the inverse method 
of separation of variables one recovers the Gaudin models as shown in refs. fs], |6|, 0. As to 
the irrational multi-parameter spectral equations, these were never discussed in the literature. 
The reason for this is that the form of the Bethe Ansatz equations determining their spectra 
drastically differs from the standard one usually obtained in the framework of r-matrix method. 
Moreover, there is no change of variables which could reduce these equations into standard form. 
This leads us to the claim that the class of completely integrable quantum systems associated 
with irrational multi-parameter spectral equations is different from the known classes of models 
and therefore its investigation is an interesting mathematical problem. 

2 The problem 

Second-order linear differential equations play an important role in many branches of mathe- 
matical physics. By an appropriate homogeneous transformation and a change of variable any 

^We mean here the so-caUed classical Yang - Baxter equations. 
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such equation can be reduced to the fohowing canonical form 



-Q^ + W{x)j^{x)=0, (2.1) 

in which W{x) and tp{x) are assumed to be analytic functions of the complex variable x. 

There are many different mathematical problems which are connected with an equation of 
the form ( |2.1| ). The simplest one can be formulated as follows: for a given function W{x) find 
the function iIj(x). The general solution of this problem obviously is 

^p{x) ^ simp ■ '4'i{x) + cos (p ■ ip2ix) (2-2) 

where i^iix) and ip2{x) are two linearly independent solutions and (p is an arbitrary mixing angle. 
The usual way to fix this angle is to impose one additional constraint on the general solution 
(|2.2| ). This constraint has to be compatible with the linearity of equation (|2.l|) and thus should 
have the form 



Comx)]=0, (2.3) 

where ^CqIV'C^)] is some appropriately chosen linear functional^. 

It is well known that, along with this trivial (and mathematically not very interesting) 
interpretation of equation , there are many others which lead to richer sets of solutions and 
are of greater theoretical importance. The essence of most of these interpretations is to allow 
some freedom in choosing the function W{x) restricting simultaneously the class of allowed 
functions ipix)- This leads to the so-called spectral versions of equation ( |2.1| ). 

Consider an example which will play a central role in our further discussion. Thereby the 
form of the function W{x) is restricted to 

N 

W{x) = Woix) + ^ enWnix), (2.4) 

n=l 

where Wo{x) and Wi{x), . . . , Wn{x) are some fixed functions and ei, . . . , e^v are arbitrary nu- 
merical parameters. Restrict the class of admissible functions ip{x) by the following -|- 1 
constraints 



Comx)] = Ciiijix)] = ... = £^[^(x)] = 0, (2.5) 

where the Cn[ip{x)], n = 0,. . . ,N are some arbitrarily chosen linear but linearly independent 
functionals. Then one can state the problem of finding those values of the parameters ei, . . . , cat 
for which equation ( |2.lD (with W{x) given by formula (|2.4D) has solutions fulfilling equations 

dU). 

It is natural to call ei, . . . , ctv spectral parameters and equation ( p.l| ) supplemented by con- 
ditions (|2.4| ) and (|2.5| ) a multi-parameter spectral equation. The set of admissible values for 
the parameters ei, . . . ,e7v we shall call the spectrum. It is easily seen that the problem (p.l|), 



(2.4), ( |2.5D is a natural generalization of an ordinary one-parameter spectral equation which 



corresponds to the case A'' = 1. 



It is not difficult to show that in general the spectrum of equations (2T), (|2.4| ), ( ^.5| ) is infinite 



and discrete. Indeed, let 'i/'i(j;, ei, . . . , e^v) and il^2{x, ei, . . . , cat) denote two linearly independent 
^In practice, one usually uses the simplest functionals: LoH'i^)] = V'(O) a-nd £oH'{x)] = ^'(0). 
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solutions of equation (|2.l|) considered as functions of the spectral parameters ei, . . . , cat. Then 
the general solution can be written in the form 



■ip{x) ~ sinc^ • ^i{x, ei, . . . , cat) + cosc^ • ^p2ix, ei, . . . , cn) (2.6) 



where ip is an arbitrary parameter (mixing angle) . Substituting ( |2.q ) into (|2.5| ) leads to a system 
of + 1 numerical equations 

lo{(p,ei,. . . ,eN) = h{(p,ei,. . . ,eN) = . . . = lN{f,ei, . . . ^cn) = 0, (2.7) 

for + 1 quantities tp and ei, . . . , cat. Here /„(c, ei, . . . , cat) denotes the value of the linear 
functional Cn[4>{x)] applied to the solution (|2.6|) . Since the number of equations coincides with 
the number of unknowns, the spectrum of equation ( |2.1| ), ( |2.4| ), ( |2.5D will be discrete in general. 
Also generally, the function ( ^.61) is transcendental which suggests that the spectrum should be 
infinite. 

It is however quite obvious that the scheme given cannot be considered a practical way for 
solving the multi-parameter spectral equations. This is so because for most functions Wq{x) 



and Wi{x), . . . , Wi^{x) the explicit form of the general solution (2^) is unknown. Thus, one can 
expect that most of the multi-parameter spectral equations of the form ( p.lj ), (p.4|), ( |2.5| ) should 
not be exactly solvable. 

Fortunately, there are exceptional cases, and these we intend to discuss in this paper. We 
think of so-called exactly solvable multi-parameter spectral equations which can be solved by 
means of purely algebraic methods and have many important applications in the theory of 
completely integrable quantum systems. One of the standard and most effective methods for 
constructing such equations is the so-called inverse method. It rests on a very simple idea: instead 
of looking for solutions of equation ( ^T| ) for given functions Wo{x) and Wi{x), . . . , Wn{x), one 
should try to reconstruct the form of these functions starting with appropriately chosen function 
V'(x). An advantage of the inverse problem in comparison with the direct one is obvious: the 
problem of solving differential equations is replaced by the problem of taking derivatives of 
known functions. However, this is only a relative simplicity: the algebro-analytic part of the 
inverse problem remains rather non-trivial which is clearly seen from the following discussion: 

Let us fix some K -\- 1 functions ipkix), k = 0, . . . , K satisfying constraints ( p.5| ) and K + 1 
sets of numbers {e^i, • • • , ekN}^ k = 0, . . . , K . Substituting these into ( |2.1D , using relation (^^), 
taking e^o = 1 and dividing finally the kth. equation by tpkix), we obtain 



dx'^ 



EefcnW'„(x) = Vfc'(x)^^|4^, k = 0,l,...,K. (2. 



n=0 

Formula (^]^) can be considered a system of K + 1 linear inhomogeneous equations for A^ + 1 
functions Wo{x) and Wi{x), . . . ,Wn{x). U K < N, then system ( |2.8D is obviously solvable. 
In this case one easily finds the explicit form of the functions Wo{x) and Wi{x), . . . ,W]\f{x) 
for which the equations ( |2.1| ), (p.^), ( |2.5| ) have K + 1 a priori known (and, hence, explicit) 
solutions. The situation changes, however, ii K > N (this is just our case, because we are 
looking for multi-parameter spectral equations having an infinite number of explicit solutions). 
In this case, the number of equations ( |2.8| ) exceeds the number of unknowns which makes the 
system ( p.^ ) overdetermined for almost all sets of functions ip^ix). The only way to get rid 
of this problem is to start with functions ipkix) for which the compatibility of the equations 



forming system (2^) would be guaranteed from the very beginning. But for this to work we 
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need some reasonable ansatz for the functions ip{x). It is hardly neccessary to emphasize that 
the problem of finding such an ansatz is far from being trivial. 

From a purely practical point of view, it is much more convenient to deal not with the 
functions iJj{x) but with their logarithmic derivatives P{x). In terms of the functions P{x), the 
ansatz which we intend to use takes an especially simple form. The substitution 

^{x) = exp I P{x')dx'^ (2.9) 



simplifies also the form of equation (|2.l| ) . The new equation 

W{x) = P^{x)+P'{x), (2.10) 
will be considered as starting point for our further considerations. 

3 Separable functions of several variables 

We shall call a function /(xi, . . . ,Xk,yi, ■ ■ ■ ,yi) oi k + I variables separable with respect to the 
variables xi, . . . ,Xk and yi, . . . ,yi if it can be represented in the form of a finite sum 

N 

/(xi, ...,Xk,yi,...,yi) = Y^ gn{xi, . . .,Xk)K{yu ■■■,yi) (3-1) 

n=l 

where hn{xi, . . . , x^), n = 1, . . . , N and gn{yi, ■ ■ ■ , yi), n = 1, . . . , N are some functions of k and 
I variables, respectively. For stressing this property we shall denote such a function by 

f{xi,...,xk,yi,...,yi) = f{xi,...,Xk\yi,...,yi) (3.2) 

Functions of several variables separable with respect to all arguments we shall call completely 
separable. Any such function can be represented in the form 

N 

,n \Xm ) (3.3) 

n=l 

where the fi^nixi), i = 1, . . . ,m are certain functions of one variable. For such functions we 
shall use the notation 

f{xi,X2, Xm) = /(2;i|x2| . . . \Xm) (3.4) 

Let us now formulate a simple lemma about separable functions. 
Lemma 1. Let r(^) be a rational function of Then 

^ + ^ = fm.) (3.5) 



and 



'^^^^ .^ =/(6ig2i6) (3.6) 



(6-6)(ei-6) (e2-C3)(e2-6) (e3-6)(e3-6) 

are symmetric and completely separable functions. 
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Proof. Any rational function can be represented as a linear combination of the so-called 
elementary rational functions having only one singularity in the complex plane. Thus, in order 
to prove separability of functions ( |3.5D and ( |3.6| ) , it is sufficient to make sure that the positions 
of the singularities of the functions /('^i|'^2) and /(^i|^2|^3) with respect to each argument do not 
depend on the values of other arguments. It is clear that the "dangerous" singularities of both 
functions (|3.5| ) and (|3.6D can only arise from their denominators. However, a simple analysis 
shows that these singularities (which, obviously, are present in each of the separate terms) cancel 
each other. The absence of these "dangerous" singularities completes the proof of the lemma. 



4 Scalar triangle equation 

Let r(^) be an arbitrary rational function of a complex variable ^. Define a complex valued 
function ^(x) by the equation 

= r(e(x)) (4.1) 
and construct from it a new function of two complex variables 

Lemma 2. The function p(x, y) obeys the following functional relations: 

pix,y)+ p{y,x) =0, (4.3) 



p{x, y)p{x, z) + p{y, z)p{y, x) + p{z, x)p{z, y) = io{x\y\z), (4.4) 
d 

-Q^p{x,y) + p'^{x,y) = uj{x\x\y), (4.5) 

with uj(x\y\z) a certain separable function. 

Proof. The proof of the anti-symmetry of the function p{x, y) immediately follows from 
definition (4.2). In order to prove the separability of the function ijo{x\y\z) we consider the 
following chain of equalities 

Auj{x\y\z) = 4[p(x, y)p{x, z) + p{y, z)p{y, x) + p{z, x)p{z, y)] = 



e{x)+e{y) e{x)+c{z) e{y)+e{z) c{y)+e{x) e{z)+e{x) e{z)+e{y) 
e(x)-^(y) e(^)-^(^) ay)-az) ay)-ax) e(^)-e(x) 



ia^) - ay)mx) - a^)) iay) - a^mix) - a^)) - amcix) - ay)) 

ria^)] I my)] ^ r[g(z)] ^^^^^ 



ia^) - ay))ia^) - az)) i^y) - az))iay) - a^)) m-amaz)-ay))' 

Lemma 1 then shows that the function u;(x|y|z) is separable. The last equation ( [4.5D can be easily 
proved if we take z = x + e m ( [4.4] ), and take the limit e — > noting that lime_^o ^pix + ^,x) = 1. 
This completes the proof. 

Hereafter we shall call equation (|4.4|) the scalar triangle equation. 



6 



5 ^-functions 

Any rational function of ^(x) and (,'{x) we shall call a ^-function. If F{x) is a ^-function, it can 
be represented in the form 

F{x)=R{ax)) + eix)G{ax)) (5.1) 

where R{£,) and G{S,) are some rational functions. The sum, difference, product and quotient of 
two ^-functions is again a ^-function. The derivative of a ,^-function is again a ^-function. We 



call a ^-function F{x) even if G(^) = in (jS.!]) and odd if = in (|5.lD . The product of 



two odd or two even ^-functions is an even ^-function, and the product of an even and an odd 
^-function is an odd .^-function. This means that the algebra of odd and even ^-functions is a 
Z2-graded algebra. The differential operator d/dx becomes in this case an odd object. 

Lemma 3. If F{x) is a ^-function, then 

a{x\y) = [F{x)-F{yMx,y) (5.2) 

is a separable function. 

Proof. Substituting ( ^.1| ) into ( ^.2[) , we can write 



where 



and 



a{x\y) = aR{x\y) + acixly) (5.3) 



^Rix\y) = ^(^^ _ g(^) {Hx)-C{y)) (5.4) 



From Lemma 1 it immediately follows that cr/j(x|y) is a separable function. In order to prove 
separability of aG{x\y), let us rewrite (|5.5[) in the form 



Ua;) - S,{y) ^(x) - ^(y) 



Using now ( [4.1] ) in the first term and applying then Lemma 1 we find that also aG{x\y) is 
separable. This proves the lemma. 

6 Bet he ansatz 

Let us look for solutions of equation ( p. 10 ) in the form 



M 



P{x) = F{x)+Y,p{x,Xi) (6.1) 



i=l 



7 



where M is an arbitrary non-negative integer, xi, . . . ,xm are still unknown numbers and F{x) 
is a ^-function. We call this form the Bethe Ansatz. Substituting ( |6.lD into ( |2.1[1| ) gives 

M 

W{x) = F\x) + F'{x) + 2 ^[F(x) - F{xi)]pix, Xi) + 

i=l 

M MM M 

+2 F{xi)p{x, Xi) + ^ p^{x, Xi) + ^ p{x, y) + Y^ p{x, Xi)p{x, Xk). (6.2) 

1=1 i=l i=l i^k 



Using formulas (4^), ( [1.5| ) and (|5.2|), we obtain 

MM M 

W{x) = F^(x) +F'{x) + 2^cj(x|xi) + Yuj{x\x\xi) + Yuj{x\xi\xk) + 

i=l i=l ij^k 

M ( M 1 

+2^p(x,Xi)< ^ p{xi,Xk) + F{xi)y . (6.3) 

We see that the first three terms in the right hand side of ( |6.3| ) represent some separable function 
of x and Xj, while the last sum of the so-called unwanted terms is, obviously, non-separable. In 
order to make the function T^(rc) separable, one should require all the unwanted terms to vanish. 
This is equivalent to the system of M equations 

M 

pixi,Xk) + Fixi) = 0, i = l,...,M (6.4) 

k=l,k^i 

which we shall call the Bethe Ansatz equations. If these equations are satisfied then 

MM M 

W{x) = F^(x) + F'{x) + 2 ^ a{x\xi) + ^ uj{x\x\xi) + ^ uj{x\xi\xk) (6.5) 

i=l 1=1 ij^k 

and hence can be written in the form 

N 

Wix) = Woix) + J2 enWnix) (6.6) 

n=l 

where Wo{x) and Wn{x), n = 1, . . . , N are certain ^-functions and e„, n = 1, . . . , are some 
numerical coefficients in which all the dependence on numbers xi,. . . ,xn is contained. We see 
that equation {2A_) takes in this case the form of a multi-parameter spectral equation. The 
role of the spectral parameters is played by the numbers ei, . . . ,eAr. The admissible values for 
these parameters are determined by the solutions of the Bethe Ansatz equations ( |6.4| ). It is 
easily seen that, for any finite M, the number of Bethe Ansatz equations coincides with the 
number of unknowns and therefore this system has a discrete set of solutions. Since this is an 
algebraic system, it has a finite number of solutions for any M, but M was an arbitrary non- 
negative integer. This means that the total number of solutions of system ( |6.4] ) is infinite and 
thus the corresponding multi-parameter spectral equation has infinite discrete and algebraically 
calculable spectrum. 



8 



7 Some simple examples 

In this section we consider three simple examples of solutions of the scalar triangle equations 
and construct the corresponding classes of .^-functions. 

Example 1. Assume that r(^) is a first-order polynomial 

r(0 = a + 6e (7.1) 

Then, from ( [4.1[) it follows that 



e'(x) = y'a + 6C(a;). (7.2) 
Solving this differential equation we obtain 

Construction of the function p{x, y) by formula ( [4.2| ) gives 

/)(x,y) = ^— . (7.4) 
x-y 



This function obeys all three relations ( [4 .31 ) - ( ^.51 ) with a trivial function uj{x\y\z): 

uj{x\y\z) = (7.5) 

In this simple case, the set of ^-functions coincides with the set of all rational functions of x. 
Example 2. Let us now assume that r(^) is a second-order polynomial: 

r(e) = a + fte + c^^- (7.6) 

Then, from ( [4.1[ ) it follows that 



i'{x) = ^a + hi{x)+ce{x). (7.7) 
Solving this differential equation we obtain 



\/ 4ac — IP' ^ /-, \ b N V 4ac — 6^ , ^/ 

= ^ smh Vc(x - t) - — , C (2;) = 77-^ — coshVc(x-t). (7.8) 

2c 2c 2^c 

Construction of function p(x,y) by formula (4.2) gives in this case 

p{x,y) = ^-cth^ix-y). (7.9) 

This function obeys all three relations ( [4. 3D - ( [4. 51 ) with a constant function uj{x\y\z): 

u;ix\y\z) = (7.10) 

In this case, the set of ^-functions coincides with the set of all hyperbolic functions of x with 
period ^-Kij^. For negative c these functions become trigonometric. 
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Example 3. Let now r(x) be a third-order polynomial of the form 

r(0 = a + 6e + + df 

Then, from (14. 11) it follows that 



Solving this differential equation we obtain 

4 c 4 

^{x) = -^V{x -t,g2,gz) - i'{x) = -V{x - t, 52, 53) 

where V{x, g2, gs) denotes the Weierstrass 7^-function with 

- 3bd 2c^ + 9bcd - 27ad^ 



92 



12 



53 



16 • 27 

Construction of the function y) by formula ( [4.2| ) gives 

= C(2; - ^,52,53) - C(a; - ^,52,53) + C(y - ^,52,53), 



(7.11) 



(7.12) 



(7.13) 



(7.14) 



(7.15) 



where C{x, g2, 93) denotes the Weierstrass ^-function. This function obeys all three relations 
( |4.3| ) - ( [4.5| ) with a non-trivial function ^(xlylz): 



w(x|y|z) = V{x,g2,g3) +'P{y,g2,g3) +'P{z,g2,g3) (7.16) 
In this case, the set of ^-functions coincides with the set of all elliptic functions of x. 

8 An equivalent description 



Note that equation ( |2.1[1| ) admits an important equivalence transformation which preserves its 
form and its spectrum. This transformation includes the change of the initial variable x 



e = ^x) 

and a linear inhomogeneous tranformation of the functions P{x) and W{x): 

m = ax) 



[£,'{x)Y 2 \di{x) 



4 \di{x) 



(8.1) 



(8.2) 



(8.3) 



In terms of the new variable ^ and new functions W and F the equation ( 2.10 ) becomes indeed 

W{i) = P\i)+P'{i). (8.4) 
If we make in (^^) the substitution 

v5'(e) 



(8.5) 
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then we obtain the transformed version of the initial hnear equation ( ^.l]) 



2 + W{0 I i^ix) = 0. (8.6) 



We see that the form of equation (|8.6| ) exactly coincides with that of the initial equation 



Linearity of the transformation (B.3) means that equation (|8.6|) is again a multi-parameter 



spectral equation having the same spectrum as (|2.1|) . The equations connected by the transfor- 
mations and ( |8.3D we shall call equivalent. 

From the examples discussed in the previous section we know that the scalar triangle equation 
has many different solutions which lead to exactly solvable multi-parameter spectral equations 
expressable in terms of rational, trigonometric, elliptic and also more complicated functions. It is 
naturally to ask, which of these equations are equivalent in the sense of the transformations ( ^.2| ) 
- ( ^.3p and which are not. For this one should solve the classification problem. The simplest way 
to do this is to find some distinguished variable ^ = (,{x) in terms of which the solutions of the 
triangle equation have a more or less unified form. The best candidate for a variable is obviously 
the function ^{x) obeying equation ( [4.11) . Indeed, in this case all functions p{x,Xi) used in the 
Bethe Ansatz ( p.l| ) transform into functions with the same denominator — where ^ = ^(x) 
and = ^(xj). Furthermore, the condition for vanishing of all unwanted (non-separable) terms 
in expression ( |6.3D transforms into the condition of regularity of this expression at the points 
^ = ^j. The most natural way to perform the neccessary calculations in this case is to start 
out immediately from the transformed version (^^) of our equation and use only very general 
informations on the form of the functions W{£,) and -P(^). This form follows from the results of 
section 5 and is given as 



WiO = AiO + ^JmBiO (8.7) 

and 



pio = «(e) + ^Jrimo (8.8) 

where A{^), B{^) and a(^), 6(^) are some rational functions and r{(,) as in ([4.11). Substituting 
(|8.7| ) and (^) into (|8.4D we obtain two independent equations for a{£^) and b{^) 



and 

BiO = b'iO + (mo + ^) KO (8.10) 



2r{0 

In the following two sections we demonstrate that these equations admit two principally different 
Bethe Ansatz solutions which we call the rational and irrational ones. 

9 The rational Bethe Ansatz solution 

Before discussing the general case, consider an important special one. Obviously, the choice 

6(0 = (9.1) 
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leads to a simpler form for system (|8.9|) , ( 8.10| ): 



A{0 = a\0 + a{0, (9-2) 

and 

B{0 = 0. (9.3) 
which is hence reduced to only one equation (|9.2| ). The Bethe Ansatz for this equation reads 

M -. 

«(0 = «(0 + E TZTT (9-4) 

i=i ^ ^» 

Substituting this ansatz into ( |9.2D and requiring regularity of the function A(^) at the points 
^j, i = 1, . . . , M we obtain 

AiO = + + 2 E ^^^^ (9.5) 

where the numbers ^j, i = 1, . . . , M have to be solutions of the system of Bethe Ansatz equations 

M , 

E F3^ + «(?0 = 0' i = l,---,M (9.6) 

k=i,k^i 



Substituting ( |9.5D and (|9.3| ) into ( ^.7| ) and using Lemma 1 we can reduce the function l^(^) to 
the form 

N 

W{0 = WoiO + E ^nWniO (9.7) 
n=l 

with Wo(0 9.nd Wn{(,), n = 1,...,A^ certain rational functions and the numbers e^, n = 
1, . . . ,N depending on the values of parameters ^i, i = 1, . . . , M which fulfill the Bethe Ansatz 
equations (|9j 



10 The irrational Bethe Ansatz solution 

The most general form for the functions a(x) and b{x) for the Bethe Ansatz for P{x) reads 

M 

«(e) = «(e) + E 7^ (10-1) 
i=i ^ ^» 

and 

^ 3- 

biO = + E TZV (10-2) 

i=i ^ ^» 

where ai,fii,(,i, i = 1, . . . , M are some unknown numerical parameters and a(^) and /3(^) are 
fixed rational functions. Substituting formulas (|10.lD and (|10.2| ) into equations (|8.9D and (|8.10| ) 
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and requiring regularity of the functions A((^) and at the points ^i, i = 1, . . . , M we obtain 

^S(^r^^""^^'^^ vm + j ^^'-'^ 

and 

where the numbers ^j, i = 1, . . . , M are again assumed to be solutions of the system of Bethe 
Ansatz equations 



M 



For the numbers ai, (3i, i = 1, . . . , M we obtain 

a^ = l, = — L=, i = l,...,M (10.6) 



Substituting (|Ta|) and into Q and using Lemma 1 we can reduce the function 

W{(,) to the form 

N 

W{0 = WoiO + E ^nWniO (10.7) 
n=l 

where Wo(0 and Wn{(,), n = 1, . . . ,N are certain irrational functions and the numbers e„, n = 
1, . . . ,N depend on the parameters ^j, i = 1, . . . , M which satisfy the Bethe Ansatz equations 

dH). 
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